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INTRODUCTION 



The brane world scenario is an impressive progress in high energy physics and cosmology in recent years. The 
' basic idea of this scenario is that the standard model particles are confined to the 3-brane, while the gravitation can 
propagate in the whole spacetime Cosmology in the brane world scenario has been widely investigated. In the 
braneworld cosmology the bulk spacetime is assumed to be Schwarzschild-Anti-de Sitter (AdS)[2] or just Minkowski 

.0- 

However, such bulks are not necessary for cosmology. In principle, we only require that the bulk space admits a 
3-dimensional maximally symmetric spacelike submanifold, which serves as our space. Besides Schwarzschild-AdS 
I^J^ (dS) and Minkowski, a 5-dimensional vacuum Taub space and the source of Taub space also admit 3-dimensional 
maximally symmetric spacelike submanifold, which are proper backgrounds for brane cosmology. We will study the 
' cosmology in the background of a sourced-Taub bulk. 
^ . The source of Taub space is a long standing problem [1] . In a recent series of works d, 0, j we successfully find 
' the source of the Taub space both in 4-dimension and higher dimensions. The components of the source region are 
, also preliminarily studied. These results strengthen the physical foundation of Taub space. 

' In the next section we briefly review the former results of the source of the Taub solution. In section III, we shall 
' explain our set up of the sourced Taub system with a moving brane. We study the cosmology of this set up in section 
IV. In section V we present our conclusion. 

o 

O ■ II. THE SOURCE OF TAUB SPACE 



' A new class of static plane symmetrj^c solution of Einstein field equation sourced by a perfect fluid was successfully 
. found in 4-dimensional spacetime in 5] . This solution is identified as the source of the Taub space since there is a 
d ' special family in this solution which can perfectly match to the Taub space. 

In 0, we generalize the 4-dimensional source of Taub solution to a higher dimensional one, which reads, 

ds^ = ~e^-^-dt^ + dz^ + e2['^-+^'^"''""'ldS2, (1) 

where 

dS^ ^ (^dx^f ^ (^^2)2 ^ ^ (rfa;"-2)2, (2) 

The above metric (jlj is an exact solution of Einstein field equation sourced by a perfect fluid, 

T^{p{z)+p{z))U®U+p{z)g„, (3) 
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where T denotes the energy momentum tensor of the fluid, U stands for 4-velocity of the fluid, g„ denotes the 
n-dimensional metric tensor, and. 



P = 



n-2 
(n-3)2 



P 



2 (n-3) 



(n - 3)2(n - 1) + 2{n - 2)26e'^^/("-3) + {n - i)fo2g2az/(„-3)l ^ 



{n - 3)(n - 1) + 2{n - 2)6e""/("-3) ^ ^2g2az/(«-3)l ^ (5) 



In 4-dimensional case, the above energy momentum can be a phantom ijj with dust and photon. The Lagrangian for 
the source is 

-Csourco = ^df^tl^d^lp -U{^p) + £dust + >Cphoton, (6) 

where the potential U{^) is rather complicated, see @ for details, £dust denotes the Lagrangian for the dust and 
'Cphoton represents the Lagrangian of the photon. The substance of n-dimensional source is almost the same as that 
of the 4-dimensional case. The reason of this similarity roots in the following fact: the only reasonable U{ip) is only 
a function of z Q . 

It is found that a family in the solution ([T]) can perfectly match to the n-dimensional Taub solution, and thus is 
the source of a n-dimensional Taub. We write the vacuum Taub metric as follows, 

ds^ ^ -z^"k^dt^ + dz^ + z^l^l^d^^, (7) 

where a — —^zji P — ;rzrr- two constants. If the vacuum region resides in z > zq and the source region 

inhabits z < zq, the matching condition at the boundary z — zq requires. 



azo 



k = ±^—, (8) 

^0 

l = ± , (9) 

and 



n — 3 

azQ = -. (10) 

n — 1 



For detailed discussions, see our original Letter Q- 



III. THE SET UP 

We consider a 3-brane imbedded in a 5-dimensional bulk. The action includes the action of the bulk and the action 
of the brane, 

S = -Sbulk + -^brane- (H) 

Here 



5buik = J d^Xy^- det (gs) (^^^5 + -Cbuikj , 



(12) 



where X = [t, z, ,x'^ ,x^) is the bulk coordinate, x^,x'^,x^ are the coordinates of the maximally symmetric space. 
M, det (55), K, i?5, >Cbuik, denote the bulk manifold, the determinant of the bulk metric, the 5-dimensional Newton 
constant, the 5-dimensional Ricci scalar, and the bulk matter Lagrangian, respectively. Note that the bulk matter 
Lagrangian £buik can be a phantom with dust and photon. 
The action of the brane can be written as, 

S'branc = / d^X^/- det(5) [kT'^ K + Lbranc) , (13) 
JM 
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where M indicates the brane manifold, Aet{g) denotes the determinant of the brane metric, Lbranc stands for the 
Lagrangian confined to the brane, and K marks the trace of the second fundamental form of the brane. x = 
(t, x^, x^, x'^) is the brane coordinate. Note that t is not identified with t if the the brane is not fixed at a position in 
the extra dimension z ^constant. We will investigate the cosmology of a moving brane along the extra dimension z 
in the bulk, and such that t is different from t. 

We set the Lagrangian confined to the brane as follows, 

-^branc = —A + L^, (14) 

where A is the brane tension and denotes the ordinary matter, such as dust and radiation, located at the brane. 
We see that our set up is a brane with tension and ordinary matter imbedded in a 5-dimensional vacuum bulk. Further 
we assume the bulk space is a sourced-Taub space. In this case, the n-dimensional metric in section II degenerates to 

55 = - dz^ - e2'(-) [{dx^f + {dx^ + [dx^f] , (15) 

where 

f{z) = az, (16) 



in the source region, and 



l{z) = az + be"^ , (17) 



22/(-) ^ k''z-\ (18) 



e^'^^) = m^z, (19) 
in the vacuum region, where a, 6, k, m are constants. Correspondingly, the matching condition ([8l [9l fTO|) 



k = ±^—, (20) 



m = ± ^ , (21) 

^0 

and 

azQ = (22) 

where the subscript denotes the value at the boundary of a quantity, and a ~ —1/2, (3 ~ 1/2. So there are only two 
free parameters, a and b, in this set-up (source and external vacuum space), just as the same of the source. This is 
natural since the source should uniquely determine the external vacuum metric. Under this matching condition, the 
metric is only at the boundary. 

By considering a brane is moving in such a bulk, we investigate the cosmology of the brane in the next section. 

IV. COSMOLOGY 

We take a method developed in We suppose that a 3-brane is moving along z direction, whose velocity is 

^ {u°,u^) ^ {i,z), (23) 

where the other 3 components of u have been omitted, since they are just 0, and a dot stands for the derivative with 
respective to the proper time of the orbit of the brane in the bulk, r. The normalization condition of u requires, 

g5{u,u)^l, (24) 
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which imphes, 



The normal of the velocity satisfies, 



and 



which yield, 



and 



e^fp -z^^l. (25) 

55(n,u) = 0, (26) 

g5(n,n) = -f, (27) 

e'^^n^i-n^i^O, (28) 



e2/(n")2 - (ni)2 = -1, (29) 
respectively. Associating (l26l) . p7)l and (124]) . we obtain 

(ni)2 = i2 + f, (30) 

(-Y = e-^^(i^ + f)J. (31) 

n can be determined up to a sign, since direction of n has not been uniquely selected. The induced metric on the 
brane reads, 

g^g5 + n(g>n = dT^ - e^'^^) [{dx^f + {dx'^f + [dx'^f] . (32) 
The second form of the brane is the Lie derivation of the induced metric along the normal direction, 

K = \cng. (33) 



whose spatial components read, 



where i, j denote the spatial index of the brane universe. By imposing a Z2 symmetry on the two sides of the brane, 
the matching condition across the brane yields. 



^2 



K^_,„ - Kg^^ = — —Sf_,^, (35) 

where lower case of Greeks denotes the index of the quantity of the brane, which runs from to 3, and s^jy represents 
the energy momentum confined to the brane, which is defined as. 



2 S{^-dct {g)L 

branc/ /o/ix 

V- det (g) 

where Lbrano is given by p4p . The spatial components of (|35p imply the Friedmann equation on the brane, 

where ptr is the density of the brane, which is defined as, 

Pbr = (38) 
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Here, the final form of the Friedmann equation is independent of the form of f{z). We note that pbr carries the 
effect of ah the matters confined to the brane, including the vacuum energy A. Only the Friedmann equation (j37p 
is not enough to determine the evolution of the brane universe. The other essential equation is the Bianchi identity 
(continuity equation), 

Pbr + iH{pbr + Pbr) = 0, (39) 

where H is the Hubble parameter, which is defined as 

H^^. (40) 
Pbr denotes the pressure of the brane for a comoving observer, 

Pbr ^s\^sl^ sl (41) 

The continuity equation (|39p also can be derived from the time component of the matching condition psp [3]. From 
(|40)) we see that the spatial coordinate -^z plays the role of the scale factor of the brane universe. The Friedmann 
equation (|37|) is a special case of the most general form of brane world model [9]. In this article, we first study a brane 
imbedded in the sourced Taub space. 

In the source region of the bulk, l{z) — az + 6e^, hence the Friedmann equation (|571) becomes, 



4z2 144a2z2(l + |e¥) 



^ = T77-T-T7p:^ir^- (42) 



In the external vacuum region, e^'^^' = m^z, the Friedmann equation (|37p becomes. 

Two comments on the Friedmann equation are listed as follows: 

1. A term l/z^, which is very alike the radiation term, but it is not the radiation, appears. This term is often 
called dark radiation in the literatures, which is the contraction of bulk Weyl tensor. In a spherically symmetric 
bulk, the physical sense of this term is the gravitational mass of the bulk space. In the present set up, it is a free 
term since z is free adjusting coordinate. This confirms the fact that the Taub space has no proper definition 
of gravitational mass since it is not asymptotically flat. 

2. The present model only admits a spatially flat universe since Taub space does not admit other type of splittings 
with a 3 dimensional space-like submanifold. It is clear that the 5-dimensional sourced- Taub (jlSp admits 6 
space-like Killing vectors, including 3 translational vectors and 3 rotational ones, which span a 6 dimensional 
Euclidean group Eq. 

We expect the present braneworld model can say something on the dark energy problem. First, we try to obtain the 
condition for cosmic acceleration in this model. Since the Friedmann equation seems fairly complicated, we deduce 
the acceleration equation by a very general function of the density of pbr and scale factor R = y/z, 

H^ = ^{pbr,R). (44) 
Associating with continuity equation (|39p . one derives 

^ = -3(p,.+p,.)^-fi?||-f20. (45) 

R Opbr oR 

One sees that the acceleration of the universe is only determined by the density and pressure iff 

C 

Cp{pbr,R)^^{pbr) + ^, (46) 



where C is a constant, otherwise the scale factor R will explicitly appear in the representation of the acceleration. It 
c 
W 



is evident that -S^ is just the spatial curvature term. Both the Firedmann equations in the source region (|42|) and in 
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the vacuum region (|43|) can not be written in the form of (l46l) , thus the acceleration has to be an exphcit function of 
R. In the source region, the acceleration becomes 



2R 



where 



2H^ 



(47) 



dR R^ \ 144a2(l + |e'^^'/2)3 



4 + 6(2 + ai?2)e''^'/^ 



(48) 



In the vacuum region, the acceleration becomes 



l = -lib(^^^'^ + ^^^'^)^^'+4^- 



(49) 



At the late-time, the dark radiation term is reasonably committed. So, if the brane is moving in the vacuum region 
of the bulk, the acceleration condition becomes 



2pbr + ipbr < 0, 



(50) 



if we require p^r > (However, for an asymptotic anti-de Sitter brane, it is incorrect since the total density may be 
negative in the late time universe. Here we do not consider this case.). It is a stronger condition than that of the 
standard case p + 3p < in the sense that it needs a smaller pressure for the same density. 
By using (fT4|) and ((38|) . we obtain 



Pbr — Pm + A, 

where pm is the density of matter confined to the brane. In the following text, we consider the case pm 
means that we consider a dust brane with tension, but without exotic matter. Substituting (|5ip into 



1 

i?4 



^ Pm 



Pn 



+ 



A2«4 



144a2i?4(i + |e'^«V2)2 72a^ Ri(l + |e'^«V2)2 I44a2^4(i + |eai?V2) 



(51) 

which 
we reach 



(52) 



However, the Friedmann equation should recovers to the standard one at the time when CMB decoupling and 
structure formation occur. ()52p cannot recover to the standard one at the early time, so it is proper to consider a 
brane moving in a vacuum bulk in the early time, where the tension of the brane is required. In the vacuum region, 
by using (|5ip the Friedmann equation (|43p becomes. 



ij2 



1 

R^ 



+ — {pl + 2p^X + X') 



Then, as usual in the brane world model, we define. 



SttG. 



eff 



288' 



(53) 



(54) 



Comparing with the previous brane models, we do not include a bulk cosmological constant since we do not find a 
proper source of Taub-(A)dS space yet. This definition of Geff is the same as the former definition in brane world 
model up to a factor Q. It has been investigated in detail except a dark radiation term, which is unimportant in the 
late time universe. 

Our scenario is that the brane moves in a vacuum bulk in the high redshift region and enters in the source region 
in some low redshift region. And the accelerating universe is driven by the source matter in the bulk. For recovering 
the standard cosmology in the middle redshift region and high redshift region, we need the brane tension. 

For convenience of numerical calculation, we write the Firedmann equation in the form of evolution with respect 
to the redshift ^. In the source region. 



jj2 

772 = -^dr(l 

Tip 



1 



(55) 
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where 



= (56) 



and 



^mp — ^jj2^ ' (57) 



- Hf' ^^^^ 



where the subscript p labels the present epoch of a physical quantity, B is defined as 



(59) 



We should not confuse the symbol of the present value p with the value of the quantity at the boundary, for which 
we use to denote. 

In the vacuum region, the Fricdmann equation keeps the same form. But B is replaced by Geff in (|54p . In contrast 
to the previous brane models, we find an interesting property of this brane world model: The equation of state (EOS) 
of the effective dark energy (defined as the ratio of pressure to energy density) can cross the phantom divide w = —1 
on a dust-brane with tension. 

The crossing —1 behavior of EOS is a deep problem and a big challenge to the fundamental physics, which is aroused 
by more accurate data: the recent analysis of the type la supernovae data indicates that the time varying dark energy 
gives a better fit than a cosmological constant, and in particular, the (EOS) w may cross —1 The dark energy 
with w < — 1 is called phantom dark energy [13], for which all energy conditions are violated. Here, it should be 
noted that the possibility that the dark energy behaves as phantom today is yet a matter in debate: the observational 
data, mainly those coming from the type la supernovae of high redshift and cosmic microwave background, may lead 
to different conclusions depending on what samples are selected, and what statistical analysis is applied |13|]. By 
contrast, other researches imply that all classes of dark energy models are comfortably allowed by those observations 



1J|. Presently all observations seem not to rule out the possibility of the existence of matter with w < —1. Even in a 
model in which the Newton constant is evolving with respect to the redshift z, the best fit w(z) crosses the phantom 
divide w = —1 jl5|. Hence the phenomenological model for phantom dark energy should be considered seriously. To 
obtain w < —1, scalar field with a negative kinetic term, may be a simplest realization. The model with phantom 
matter has been investigated extensively [1^], and a test of such matters in solar system, see However, the EOS 
of phantom field is always less than —1 and can not cross —1. It is easy to understand that if we put two scalar 
fields into the model, one is an ordinary scalar and the other is a phantom: they dominate the universe by turns, 
under this situation the effective EOS can cross —1 It is worthy to point out that there exist some interacting 
models, in which the effective EOS of dark energy crosses —1 [l9||. More recently it has been found that crossing — 1 
within one scalar field model is possible, the cost is that the action contains higher derivative terms [2fl|] (see also 
[2l|). Also it is found that such a crossing can be realized without introducing ordinary scalar or phantom component 
in a Gauss-Bonnet brane world with induced gravity, where a four dimensional curvature scalar on the brane and a 
five dimensional Gauss-Bonnet term in the bulk are present (2^ . 

To explain the observed evolving EOS of efi^ective dark energy, we calculate the equation of state w for the effective 
"dark energy" caused by the tension and term representing brane world effect by comparing the modified Friedmann 
equation in the brane world scenario and the standard Friedmann equation in general relativity, because all observed 
features of dark energy are "derived" in general relativity. Note that the standard Friedmann equation in a four 
dimensional spatially fiat universe can be written as 

= —{Pm+Pde), (60) 

where the first term in RHS of the above equation represents the dust matter and the second term stands for the 
effective dark energy, and G is the Newton constant. Comparing ((60| with (|37p . one obtains the density of effective 
dark energy, 

Pde^--r^ + — rrw: ^p™- 61 
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0.2 0.4 0.6 0.8 1 1.2 1.4 



FIG. 1: The evolution of the density for the effective dark energy p^e with respect to the redshift ^, where pcri denotes the 
present critical density. The left part (low redshift region) is the source region, and the right part (high redshift region) is the 
vacuum region. In this figure a — —1.2. The boundary of the source region inhabits at Ro — 0.645. The brane sweeps across 
the boundary when ^ = 0.549. 



Since the dust matter obeys the continuity equation and the Bianchi identity keeps valid, dark energy itself satisfies 
the continuity equation 

^+3i/(p,e+Pe//)=0, (62) 

where Pef / denotes the effective pressure of the dark energy. And then we can express the equation of state for the 
dark energy as 

Peff -.1 dlnpde 
Pde 3dln(l+0 

Clearly, if is greater than 0, dark energy evolves as quintessence; if ^"^^"j^^'^ is less than 0, it evolves as 

phantom; if fi^^^i^'^-^ equals 0, it is just cosmological constant. In a more intuitive way, if pde decreases and then 

increases with respect to redshift (or time), or increases and then decreases, which implies that EOS of dark energy 
crosses phantom divide. The more important reason why we use the density to describe property of dark energy 
is that the density is more closely related to observables, hence is more tightly constrained for the same number of 
redshift bins used ^24j. 

We find a concrete example in which the EOS of the effective dark energy crosses —1. For the numerical improve- 
ment, we should nail down the parameters in this model, which satisfy the requirements of theory and observation. 
For the theoretical aspect, the boundary condition of the source region yields 



QUU 

k 



azQ 



b^-Ae"—, (64) 

and ([22|) . We only consider the large branch in 0, hence h is fixed by ([22| and ((64|l . Therefore, to this brane model 
imbedded in sourced Taub background itself, we have only one free parameter. It is a or zq, which implies the thickness 
of the source region. The partition of dark radiation should be smaller than that of cosmic microwave background 
(CMB), otherwise the universe will bounce at some high redshift, thus Q.dr < 10^''. Here we set fidr = 10^^, Rp = 1, 
^mp — 0.27 and U\ — 1.408. The present part of baryon density is approximately 4% [2^, and we have several 
evidences of the existence of the non-baryon dark matter, which are independent of ACDM model [2^. Here we take 
ilmp as the same value of the best fit of ACDM model, fix is derived from ([55]) at ^ = 0. Figure 1 illustrates the 
evolution of the effective density. It is clear that the EOS of the dark energy crosses — 1 at ^ = 0.2 ~ 0.3. This point is 
also confirmed by the direct plot of the EOS of effective dark energy in figure 2, which well explains the observations 

The most significant parameters from the viewpoint of observations is the Hubble parameter H{z), which carries 
the total effects of cosmic fluids. Except the indirect data of H{z), such as the luminosity distances of supernovae, 
the direct H(z) data appear in recent years, which can be used to explore the fine structures of the Hubble expansion 
history [27j . There is an important new feature of H{z) data which is not implied by the previous indirect observations 



9 




0.1 0.2 0.3 0.4 0.5 



FIG. 2: The evolution of Wde with respect to ^. The point of crossing —1 (C-point) appears at ^ = 0.2 ^ 0.3. The parameters 
are the same as Figure 1. 
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FIG. 3: The evolution of the Hubble parameter with respect to redshift. The unit of H is km s ^ Mpc ^ . The parameters are 
the same as Figure 1. 



of Hubble parameter: It decreases with respect to the redshift ^ at redshift ^ ~ 0.15 and ^ ~ 1.5, which means that 
the total fluid in the universe behaves as phantom. We find that the present brane model also can partly realize this 
property with the same parameters of the example in which the EOS of the effective dark energy crosses —1. Figure 
3 illustrates the evolution of the Hubble parameter. At ^ ~ 0.15, there is a clear drop on H{z) diagram in Figure 3. 

Figure 1 and Figure 3 illuminate that there are cusps at the boundary of vacuum and source region, which is related 
to the fact that the metric on the boundary is only C^, while second derivatives with respective to coordinates are 
involved in Einstein equation. 



V. CONCLUSION 



We investigate the brane world model in the sourced-Taub background. The Friedmann equations are obtained 
both in the source region and the vacuum region. 

The recent observations imply that the EOS of dark energy may cross —1 and there may be some drops on the 
Hubble diagram. These two phenomena are serious challenges to the physical cosmology. In the braneworld model 
with a sourced Taub background, for a reasonable parameter set the EOS of the effective dark energy can cross — 1 
in when the brane is moving in the source region. And furthermore, the Hubble parameter has a drop in the low 
redshift region (the source region). 

In the present work, we do not include a bulk cosmological constant and thus the brane tension can not be too 
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large. When the brane is moving in the vacuum region, the quadratic term of p will dominate the linear term p before 
long. In the high redshift region of the vacuum bulk, the model is treated as a toy model. Our main aim of this 
article is to study the behavior in the source region. The detailed evolution of the brane universe in sourced Taub 
background in the high redshift region needs to be investigated further in the future. 
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ported by grant No. C00022 from the Korea Research Foundation. 
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